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ABSTRACT 

The helium-merger gamma-ray burst progenitor is produced by the rapid accretion onto a com- 
pact remnant (neutron star or black hole) when it undergoes a common envelope inspiral with its 
companion's helium core. This merger phase produces a very distinct environment around these out- 
bursts and recent observations suggest that, in some cases, we are detecting the signatures of the past 
merger in the GRB afterglow. These observations allow us, for the first time, to study the specific 
features of the helium merger progenitor. In this paper, we couple population synthesis calculations 
to our current understanding of gamma-ray burst engines and common envelope evolution to make 
observational predictions for the helium-merger gamma-ray burst population. Many mergers do not 
produce GRB outbursts and we discuss the implications of these mergers with the broader population 
of astrophysical transients. 

Subject headings: Supernovae: General, Stars: Neutron 



1. INTRODUCTION 

The black hole accretion disk gamma-ray burst (GRB) 
engine has be come standard for lon g- and short-duration 
GRBs alike (|Popham et al.l 119990 . For long-duration 
bursts, the most commonly-invoked progenitor is the 
collapse of a rapidly rotating massive star down to a 
black hole. But achieving the high angular momentum 
needed to produce a disk around the newly formed black 
hole has proved problematic and a number of progen- 
itors have been inv oked to provide these high angular 
mome n tum profiles dFrver fc Wooslevl 119981; IFrver et all 
19993 lYoon fe Langerl 120051: iWooslev fe Hegerl 120061: 
Belczynski et al. 20071: llvanova fc Podsiadlowskil 



20031: iPorteeies-Zwart. Dewi. fe Maccarom 
van den Heuvel fc Y oon 120071 : IfVveT et all 
Podsiadlowski et al.ll2010h . 

One such model invokes the merger of a compact 
remnant with its companion. As the compact rem- 
nant spirals into the center of the helium core of 
an evolved companion, the accretion rate can exceed 
0.01 M f^s -1 , enough to pow er a strong gamma-ray 
burst (jFrver fc Wooslevl 119981 ). The basic idea behind 
this model is that as the compact remnant spirals into 
its companion, the orbital energy lost both spins up the 
helium core and ejects the hydrogen envelope. If the 
compact remnant is a neutron star, the rapid accretion 
in the inspiral is likely to cause it to collapse to a black 
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hole. In this manner, the helium merger model provides 
a natural way to ensure enough angular momentum to 
form a disk around a central black hole. 

Simulations by Zhang & Fryer (2001) showed that 
the inspiral would produce both the high angular mo- 
menta and accretion rates needed to create a black 
hole accretion disk engine. Although there is some 
concern that helium merg ers produce disks wit h too 
much angular momentum ( DiMatteo et al.l [2002',) , cur- 
rent analyses find this angu lar momentum is appropriate 
for accretion disk engines ( Komissa rov fc Barkovl 120091 : 



mgi: 

iBarkov fc~K omissa rovll2011 ). The helium- merger model 
predicts a feature that may be observable in the resultant 
explosion: the existence of a shell (or torus) of merger- 
ejecta surr ounding the burst . 

Initially, IFrver fc Wooslevl (|1998l ) believed the helium 
merger model would be difficult to distinguish from nor- 
mal GRBs because the ejecta from the common enve- 
lope phase would lie along the orbital plane (the same 
plane as the disk) and there would be very little inter- 
action between the ejecta and the GRB jet. This pic- 
ture changed with the obs ervation of GRB12 1225, the so- 
called "Christmas Burst" (|Thone et al.ll2^TTI ). This burst 
exhibited blackbody emission arguing for a compact pro- 
genitor (~ 10 11 cm) followed by a second blackbody com- 
ponent with characteristic radius of ~ 10 14 cm. Simula- 
tions of the jet interactions with more realistic common 
envelope ejecta profiles (also based on recent simulations) 
found that the helium merger model could produ ce the 
thermal emission in this burst (| Thone et al.l 12011). The 
X-ray Flare 060218 (and perhaps a few other low-redshift 
bursts) may have similar characteristics, suggesting that 
these helium merger events may be identifiable after all. 

To truly compare these systems to the observed data, 
we must understand the entire population of compact 
merger events. In this paper, we undertake a series of 
population synthesis models studying the characteristics 
of compact merger events. Section [2] describes the pop- 
ulation synthesis calculations and presents the basic re- 
sults of these calculations: rates and characteristics of 
merging systems. In section [31 we analyze these results, 
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making predictions for luminosities and ejecta positions 
for merger outbursts. We conclude with a discussion of 
the potential observational outcomes of this black hole 
accretion disk progenitor. 

2. NEUTRON STAR MERGER POPULATIONS 

2.1. Population Synthesis Calculations 

To study the helium merger progenitor, we use 
the population synthesis code StarTrack to calcu- 
late the numbers and properties of merging sys- 
tems. The full description of the co d e ca n 
Belczvnski. Kalogera. fc Bulikl (|2002D ; 
(12008ft. The code utilizes a set of stel- 



be found in 



iBelczvnski et al 



lar models (|Hurlev. Pols, fc Toutl (|2000D ; slightly modi- 
fied from its original version) that allow for the evolution 
of stars at different metallicities. The model for compact 
object formatio n adopted in the co de has been signifi- 
cantly revised ([Dominik et alj 12012). During core col- 
laps e, fallback and direct BH formation is now accounted 
for (| Fryer fc Kalogera 2001) and the newly born objects 
receive natal kicks (Hobbs et al. 2005). The formation of 
low mass NSs throu gh electron capture s u perno vae is also 
accounted for (e.g., Podsia dlowsk et all (|2004| )). Binary 
interactions are treated in detail, and the various physical 
processes have been calibrated using either results of de- 
tailed evolutionary calculations (e.g- JWellstein fc Langerl 
(fl999h for mass transfer sequences), or specific set s of ob - 
servations (e.g., iLevine. Rappaport. fc Zoicheskil (|2000ft 
for tidal interactions). The mass loss prescription for 
winds in massive stars, the stellar binding energy used 
in common envelope evolution, and the compact rem- 
nant mass formation pres cription have all been up- 
dated (|Dominik et al.ll2012f ). 

2.2. Rates 

In each calculation, we have evolved N = 2 x 10 5 bina- 
ries with the primary mass range Mi = 5 — 150 M and 
secondary mass range Mi = 0.08 — 150 M Q . The stars 
were evolved with solar metallicity (Z — 0.02). Each bi- 
nary system is initiated by four parameters, which are as- 
sumed to be independent: the primary mass M\ (the ini- 
tially more massive component), the mass ratio q = j^, 
where M2 is the mass of the secondary, the semi-major 
axis a of the orbit, and the orbital eccentricity e. 

For both, single stars and binary system primaries, we 
use the initial mass function adopted from Kroupa & 
Weidner (2003), 
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0.08 < M < 0.5 M 
0.5 < M < 1.0 M© 
1.0 < M < 150 M 



(1) 



where parameter a; m f = 2.35 — 3.2, with our standard 
choice of 2.7. 

We adopt a flat mass ratio distribution that is con- 
sistent with the recent observational results (Kobulnicki, 
Fryer & Kiminki 2006), 



= 1 



(2) 



in the range q = — 1. Given value of the primary mass 
and the mass ratio, we obtain the mass of the secondary 
M 2 = qM\. 



TABLE 1 
Population Synthesis Models 



Model 


Ace 


Hertzprung Gap 


Metallicity 






Fatc a 




Modell 


physical A b 


CE formulae 


solar 


Model2 


3 X physical A 


CE forumulac 


solar 


Model3 


constant A = 1 


CE formulae 


solar 


Model4 


physical A 


HG always merges 


solar 


Model5 


physical A 


CE formulae 


10% solar 


Model6 


physical A 


HG always merges 


10%solar 



a Stars in the Hertzprung gap do not have strong entropy 
gradients dividing the helium core and hydrogen envelope. 
Without these sharp gradients, every CE phase will proceed 
through the merger with the helium core. We choose two 
fates: ignore this entropy argument ("CE formulae") or in- 
clude this effect and assume on Hertzprung gap CE scenarios 
lead to a merged system ( "HG always merges"). For more 
details, see IDomini k et al. (20l3). 

Physical A c orresp onds to the IXu fc Lil pTiTJ ) used by 
IDominik etaO 120121) . 

The distribution of the initial binary separations is as- 
sumed to be flat in the logarithm (Abt 1983), 
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where a ranges from a minimum value, such that the pri- 
mary fills at maximum 50% of its Roche lobe at ZAMS, 
up to 1O 5 R . 

We adopt the thermal-equilibrium eccentricity distri- 
bution for initial binaries, 



H(e) = 2e, 



(4) 



in the range e = — 1 (e.g., Heggie 1975; Duquennoy & 
Mayor 1991). 

We have evolved only binaries massive enough to po- 
tentially produce at least one compact object and there- 
fore produce a helium merger. If we extend the evolved 
population to include the entire mass spectrum and if 
we include single stars as well (with binary fraction 
/bi = 0.5), each calculation with a lm { = —2.7 would 
correspond to a stellar mass formed in a starburst mass 
of 4.5 x 10 7 M . 

Table[2]shows the total number of he-mergers produced 
using 2 x 10 5 binaries. The rates are divided into 3 cat- 
egories: main-sequence merger and mergers with small, 
and large, helium core masses. These rates correspond 
to the rate per 4.5 x 10 7 M of star forming mass. The 
number in parantheses show the number of the total sys- 
tems produced because of neutron star kicks. Although 
kicks produce many of the merging systems (30%), kick- 
induced mergers account for a smaller fraction (~10%) 
of low-mass helium core mergers and almost nobe (<1%) 
of the massive helium mergers believed to be GRBs. 

These rates can be folded with any desired star forma- 
tion history to give actual time dependent helium merger 
rate in a given galaxy. Table [3] shows the merger rate per 
Myr in a Milky- Way massed galaxy as well as the vol- 
umetric (Gpc -3 yr -1 ) rate. If we wish to compare to 
GRBs, a rough estimate of the rates should include a 
combination of solar and low metallicity rates (e.g. 50% 
each of Model 1 and Model 5 or 50% each of Model 4 
and Model 6). The resultant rates predicted are 629, 857 
Gpc -3 yr -1 for respective models assuming normal CE 
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TABLE 2 
Merger Rates (per 4.5 x 10 7 M©) 



Model He-Merger He-Merger Rate Main-Sequence 

(M He < 4.0M©) (M Ke > 4.0M©) Merger Rate 



Modell 
Model2 
Model3 
Model4 
Model5 
Model6 



718(110) a 
507(35) 
967(94) 
929(24) 
1505(1) 
1790(1) 



939(0) 
406(0) 
227(3) 
1020(2) 
2019(9) 
3033(6) 



579(175) 
1267(411) 
1566(478) 
580(175) 
796(249) 
755(262) 



a The left number in each column denotes the total number of 
systems produced in each category. The numbers in parantheses 
show the number of these systems formed due to neutron star 
kicks. 



TABLE 3 

Merger Rates II (Myr -1 , 



Gpc" 



- 1 ) 



Model 


He-Merger 


He-Merger Rate 


Main-Sequence 




(M Ho < 4.0M©) 


(M Ho > 4.0M©) 


Merger Rate 


Modell 


30.8(308.2) a 


40.3(403.0) 


24.9(248.5) 


Model2 


22.0(220.3) 


17.6(176.4) 


55.0(550.5) 


Model3 


43.2(431.5) 


10.1(101.3) 


69.9(69.9) 


Model4 


39.9(399.2) 


43.8(438.3) 


24.9(249.3) 


Model5 


63.7(636.7) 


85.4(854.2) 


33.7(336.8) 


Model6 


75.2(752.4) 


127.5(1274.8) 


31.7(317.3) 




20 40 60 HO 100 

Stellar Companion Mass (M e ) 




Stellar Companion Mass (M e ) 



a The left number in each column denotes a rate per Milky- Way 
massed galaxy (Myr -1 ). The numbers in parantheses show vol- 
umetric rate (Gpc -3 yr _1 ). 

formulae for Hcrtzsprung gap mergers versus immediate 
mergers for these Hertzsprung gap mergers. 

The Fryer et al. (1999) Galactic rate for their standard 
model is 35 Myr" 1 , very close to our rate of high-mass 
helium cores. Such high rates mean that even if only 
10% of these systems produce GRB jets, they can ex- 
plain a large fraction of the GRB rates. The low-mass 
helium star merger and main-sequence merger rate is not 
significantly higher and, assuming these mergers produce 
peculiar supernovae, they can account for no more than 
1% of all core-collapse sueprnovae. 

Like long-duration GRBs, the delay time from star for- 
mation to merged event is limited to the stellar evolu- 
tion time (in the he-merger case, it is the stellar evolu- 
tion time of the companion star). For massive helium 
mergers, the merger occurs within 20 Myr of the initial 
star formation time (for main-sequence star mergers, the 
distribution can extend to 100 Myr). These objects will 
likely trace star formation. 

2.3. Mass and Orbital Separation Distributions 

Our populations demonstrate that compact remnant 
mergers can occur for a wide range of star characteristics: 
mass, evolutionary phase and, hence, orbital separation 
at the time of the merger. In section [3l we will use these 
characteristics to determine observable features of these 
mergers. Before we do so, let's review the characteristics 
of these populations. 

Figures ll!2l show the distribution of companion masses 
for 4 different merger populations: those with hydrogen 
cores merging through mass transfer, those with hydro- 
gen cores merging when the compact remnant is kicked 
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Fig. 1. — Companion star mass distribution for 4 different popu- 
lations: those with hydrogen cores merging through mass transfer, 
those with hydrogen cores merging when the compact remnant 
is kicked into the companion, those with helium cores merging 
through mass transfer, and those with helium cores merging when 
the compact remnant is kicked into the companion. Except for 
helium merger systems (which have companion masses ranging up 
to 60 M©, often with a secondary peak at these high masses), most 
of the systems have companion masses peaking strongly at low 
masses. 
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Fig. 2. — Companion star mass distribution for 4 different pop- 
ulations (see fig. [T}. Note that for low metallicity systems, the 
helium merger population has a much broader high-mass tail, pro- 
ducing mergers with helium cores in excess of 90 M©. 

into the companion, those with helium cores merging 
through mass transfer, and those with helium cores merg- 
ing when the compact remnant is kicked into the com- 
panion. As we might expect with our steep initial mass 
function, the companion masses are peaked to lower mass 
(below 20 M©). However, the distribution is broad, ex- 
tending beyond 60 M Q . He-merging systems have a sec- 
ondary peak at higher masses and the most massive he- 
lium cores will come from this peak. For low-metallicity 
stars, this secondary peak occurs nearly at 90 M . 

Similarly, the distribution of companion radii is ex- 
tremely broad (Fig. [3]). As expected, the evolved stars 
(with helium cores) can be much larger (up beyond 5 
A.U.). The helium core sizes tend to lie in a narrow 
range between 0.5-3 R Q . 

The wide variation in these systems will produce a 
wide variation in the observational features of these ex- 
plosions. 

3. OBSERVATIONAL COMPARISONS 

The characteristics of the compact binaries (compan- 
ion masses, evolutionary phases, radii and orbital sep- 
arations) allow us to estimate both the explosion lumi- 
nosities and spectral features. 

3.1. Luminosities 

For explosion luminosities, we will divide the discussion 
into two categories based on the accretion rate. If the sys- 
tem has a massive helium core (> 4M Q ), the accretion 
is sufficient to power a strong black-hole accretion disk 
engine discussed in classical GRBs. In these cases, even 
the neutrin-annihilation powered engine may be strong. 
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Fig. 3. — Companion star radius distribution for 4 different pop- 
ulations in Figure [T] The evolved stars (with helium cores) have 
much higher radii (up above 5 A.U.). Because there is very little 
variation between models, we show only 4 of our models. 
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These high accretion-rate systems fit the classic Helium 
Merger model proposed by Fryer & Woosley (1998). We 
then study low-mass helium core and hydrogen core sys- 
tems. 

3.1.1. Massive Helium-core Merger: GRB-like 

Our population synthesis models predict a wide set of 
merging systems involving a compact remnant (neutron 
star or black hole) with its companion. As we showed 
in section 12.21 8-50are mergers with evolved stars with 
helium core masses above 4M . With a rate of 10- 
120 Myr -1 , mergers with massive helium cores can make 
up a significant fraction of our GRB population. In this 
section, we use our current understanding of disk-engine 
models to predict explosion luminosities from the helium 
core mass distribution. 




Fig. 4. — Mass distribution of helium cores in our population of 
he-mergers. The two largest peaks exist for low-mass helium cores: 
1.4 and 3.4 Mq. But the distribution is widespread (half of all 
merging systems involve helium cores with masses above 4 Mq ) . 

Zhang & Fryer (2001) studied the accretion onto a 
neutron star during an inspiral. They found that, the 
more massive the helium core, the higher the accretion 
rate (and, ultimately, the more energetic GRB explo- 
sion). Our population synthesis models provide us with 
the distributions of helium core masses at the onset of 
the merger is shown in Figure HI Although the helium 
core masses peak at low masses, the distribution of he- 
lium core masses for solar metallicy models extends be- 
yond 25 Mq and roughly half of the systems have helium 
cores with masses above 4M©. At low metallicities, the 
maximum helium core mass extends from 25 M up to 
45 M . 

We can estimate the accretion rate onto the compact 
remnant (M rcm ): 

M Icm = ^AiHcAfrcc-fall = 

77[2Af Hc \/2.0G(M He + M rcm )}/[nR Hc ] (5) 

where if rC e-faii is the free-fall timescale, Mh c , ^?He are 
the helium core mass and radius respectively, M Tem is 
the compact remnant mass and G is the gravitational 



constant. The accretion rate is less than the typical free- 
fall time because the material has thermal pressure and 
considerable angular momentum. To account for this, we 
include an efficiency parameter r]. Fitting to the Zhang & 
Fryer (2001) results, we find the value of r\ to be roughly 
15%. 




0.001 - 7 

10 

He-core Mass (M ) 

Fig. 5. — Average accretion rate as a function of helium core 
mass. The dots show results using the helium cores from the 
Woosley et al. (2002) models. The accretion rate data is fit by 
a 2-piece power law of the helium core mass. Below 4Mq, the 
power increases dramatically, driving the accretion rate down sig- 
nificantly as the helium core mass decreases. 

By using the helium core radii as a function of mass 
from the Woosley et al. (2002) pre-collapse helium cores, 
we derive an accretion rate as a function of helium core 
mass is shown in figure [5] These helium cores are evolved 
(at collapse), but their radii are not too different from 
the unevolved cores used by Zhang & Fryer (2001). This 
mass loss rate is fairly well fit by a 2-part power law: 

log(M rcm ) = 2 log(M H c) - 3.75 if M Hc > 4.0 

5 log(MHc) — 5.55 otherwise. (6) 

For helium cores with masses below 4 M@ , the accretion 
rate drops rapidly. These low-mass cores are less likely to 
drive strong explosions, and for the purposes of our dis- 
cussion of helium mergers, we will focus on those helium 
cores with masses above 4 Mq . 

Zhang & Fryer (2001) estimated the explosion lu- 
minosity as a function of accretion rate for a he- 
lium m erger GRB assuming b oth the Blandford- 
Znajek (jBlandford fc Znajeld Il977t ) and neutrino anni- 
hilation mechanisms. For neutrino annihilation, we use 
the formula from Zhang & Fryer(2001): 

log[ J L t/ ,P(ergs- 1 )]w 43.6 + 4.89 loglAf^/^.OlMos- 1 )] 
+3.4a rom (7) 

where a rcm = (J rem c)/(GM r 2 em ), J rcm is the angular mo- 
mentum of the remnant and c is the speed of light. Val- 
ues of a rom for these mergers tend to lie in the range of 
0.75-0.95. We will assume 0.9 for our estimates. This 
estimate assumes a roughly 3M Q black hole which, for 



6 




He-core Mass (M ) 



Model 1 

L Model 2 

> ., - Model 3 

| ! Model 4 

Model 5 

Model 6 



J'l - - - r l 

! L i.J ; 



-I |l 



46 48 50 

Log Luminosity (erg s -1 ) 



Fig. 6. — The luminosity versus helium core mass for our two 
power sources: neutrino annihilation (L v ,n), Blandford-Znajek 
emission (Lbz)- The stronger dependence of neutrino annihila- 
tion on the density and temperature (and hence accretion rate) 
means that the luminosity from neutrino annihilation is extremely 
sensitive to the helium core mass. 

the purposes of this study, is appropriate. However, 
other formulae exist that will produce different energy 
distributions, e.g. Zalamea & Beloborodov (2011). Our 
prescription produces a broader energy distribution than 
Zalamea & Beloborodov (2011) and, to demonstrate the 
range of energy distributions, we will focus on these for- 
mulae. However, one should bear in mind that the distri- 
bution of neutrino annihilation outbursts may be more 
narrow than predicted by our results. 

The comparable equation for Blandford-Znajek emis- 
sion is (Pop ham et al.lll999f ): 

log[ J L B z(ergs- 1 )] « 50.0 + 2 log[a rem M renl B Blag ] (8) 

where B mag is the magnetic field strength. As with 
Zhang & Fryer (2001), we will assume the magnetic en- 
ergy density is equal to the kinetic energy density: 

B 2 = Anpv 2 a M Icm c/r 2 g (9) 

where we have approximated the density (p) and velocity 
(v) of the material using the values at the event horizon: 
r g = 2GM rcm /c 2 . Again, we set a rom to 0.9 for our 
estimates. The corresponding luminosity for both these 
models as a function of helium core mass is shown in 
figure [51 Neutrino annihilation is much more sensitive to 
the helium core mass and leads to a broader spread of 
explosion energies. 

These formulae allow us to estimate the luminosity 
of the outburst from our helium merger population. 
The distributions for both the neutrino annihilation and 
Blandford-Znajek models are shown in figure [7] As 
expected, the sensitivity of the neutrino annihilation 
emission model leads to a broad spectrum of explosion 
luminosities. Our simplified Blandford Znajek model 
predicts a relatively narrow range of luminosities from 
10 49 - lO^ergss- 1 . 

The merger process ejects the hydrogen atmosphere of 
the star in a shell of material with an density enhance- 
ment along the orbital plane. It was this shell that Thone 
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Fig. 7. — The distribution of luminosities from the helium merg- 
ers from our population synthesis calculations. The neutrino an- 
nihilation model's strong sensitivity on the helium core mass leads 
to a broad distribution in the explosion energies. 

et al. (2011) exploited to explain the peculiar features 
of the Christmas burst. They argued that interactions 
of the jet with this shell produced the late-time, large 
radius blackbody emission. We estimate the outer ex- 
tent of this shell using the same analysis of Thone et al. 
(2011): assume the inspiral occurs in ~3 orbits and the 
material ejected in this common envelope phase moves 
outward at the escape velocity. Although these two as- 
sumptions are not strictly true, we can not yet model 
this process exactly (e.g. Passey et al. 2011, Ivanova et 
al. 2012). These assumptions approximate the ultimate 
result within a factor of 2 in most cases. Using these 
assumptions, the corresponding distribution of shell po- 
sitions for our populations are shown in figure [5] 

Our population helium mergers predicts a range of lu- 
minosities versus maximum shell radii (Fig. There is 
considerable scatter in the relationship between the outer 
shell radii and luminosity. But there is a trend that the 
most luminous explosions have the furthest shell distri- 
butions. 
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Fig. 8. — The distribution of outer positions of the helium merger 
ejecta from our population synthesis models. 

3.1.2. Supernova-like Explosions: Low-Mass Helium Cores 
and Main Sequence Mergers 

If the compact object merges with a lower-mass helium 
core or a main-sequence star, the accretion rate onto the 
remnant will be much lower. For these systems, the black 
hole accretion disk engine is unlikely to be valid and it 
is more difficult to produce even a rough estimate of the 
explosion energy. 

Before the common envelope phase, many of our sys- 
tems are likely to be X-ray binaries. At the onset of the 
common envelope phase, the X-ray binary can evolve, 
becoming more powerful. The system may even begin 
to produce strong jets. These so-called, micro-quasars, 
have been observed in numerous systems (see Mirabel 
2007 for a review). 

What happens as the compact remnant spirals into its 
companion is more difficult to determine. In these merg- 
ers, the accretion rate is slow enough that the neutron 
star may not initially collapse to form a black hole. The 
accretion rate remains high enough that any neutron star 
magnetic field is likely to be buried. Even without mag- 
netic fields, these accretion even ts can produce explosions 
on pa r with normal supernova (jFrver et alJl2006at iFrverl 

\mm . 

Presumably, once the compact remnant collapses to 
a black hole, the accretion disk will produce a jet. In 
these hydrogen or low-mass helium cores, the accretion 
rate is insufficient for the neutrino annihilation model 
to work, but we still do not understand the black hole 
accretion disk mechanism well enough to predict whether 
such systems will produce jets with explosions at GRB 
levels or at sub-supernova levels. 

3.2. Nucleo •synthetic Yields 

Two sites have been identified in long-duration GRBs 
to produce 56 Ni: explosive nucleosynthesis in the shock 
(this shock burning is what produces the 56 Ni in su- 



pernovae), or disk- wind nucleosynthesis (|Surman et al.l 

2oo6i [Iml). 

Explosive nucleosynthesis produces 56 Ni in the silicon 
and oxygen layers of an evolved core in normal super- 
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Fig. 9. — Luminosities versus outer shell radii for our helium 
merger population. Our Blandford-Znajek, magnetic field, engine 
(crosses) tends to produce stronger explosions in general, but the 
strongest explosion energies arise from our neutrino annihilation 
model (circles). The strongest energies tend to have the widest 
separations, but a there is considerable scatter in the maximum 
radius/luminosity relation. 
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Fig. 10. — The distribution of outer positions of our low-mass 
helium cores and hydrogen stars. The distribution is much broader 
than the high-mass helium cores, with many more systems with 
quite low outer shell radii. The broad range of shell separations 
from these systems is in large part due to the fact that a sizable 
fraction of these systems (30% of hydrogen stars, 17% of low-mass 
helium stars) arise from systems where the compact remnant is 
kicked into its companion. 



novae or classic collapsars. Without an evolved core, 
it is difficult for the helium merger model to produce 
significant 56 Ni from shock interactions. Unless con- 
siderable material is processed and ejected in the disk, 
very little 5 6 Ni is produced. Like the f allback GRBs 
discussed in iFrver. Hungcrford fc Young! ()2007bl ) , these 
GRBs could produce low-nickel yield outbursts. Al- 
though this may not affect the peak of the light-curve (in 
most of these systems, shock heating dominates the light- 
curve at peak: Fryer et al. 2007b, Fryer et al. 2009), it 
does mean that if explosive nucleosynthesis is the domi- 
nant source of 56 Ni, these explosions will not have typical 
56 Ni decay tails at late times. 

However, if disk winds dominate the 56 Ni production, 
the nucleosynthetic yields of GRB-like he-mergers will be 
indistinguishable from collapsar long-duration GRBs. 



3.3. Spectra 

Contamination from the ejected hydrogen envelope is 
the strongest constraint for he-mergers. However, the 
hydrogen distribution in the helium merger model than 
in typical stellar envelopes. Even though the hydrogen 
remains within f0 15 cm of th e central engine, t he dis- 
tribution is very asymmetric. iPassv et al.l (|2012t ) found 
that 90% of the mass ejected was distributed within 30° 
of the equator. It is more difficult to accurately describe 
the distribution of the remaining 10%. Combined with 
the fact that the explosion is also asymmetric, it is pos- 
sible to sweep up very little mass. 

How low must the swept-up hydrogen ma ss be? For 
symmetric explosions. iHachinger et alT (j2012f ) found that 
roughly 0.1 solar masses of hydrogen could be hidden in 
the explosion. For asymmetric supernovae, this number 
may be a bit greater. This would require only 1-2% of the 
stellar-envelope hydrogen remain within the supernova 
opening angle (perhaps within 30° of the rotation axis). 

To determine how well the low density funnel can fo- 
cus th e explosion we have simulated usi ng th e SNSPH 
code (jFrver et al.l l2006d : lEllinger et all l20ll a 5.6 x 
10 51 erg explosion focused in a 30° wedge of a 25 M 
where the density in this wedge is lowered by a factor of 
100. The explosion at 0.47d is shown in figure[TTJ In this 
calculation, roughly 0.0001 M Q of hydrogen is swept up 
in the explosion. This hydrogen mass is 3 orders of mag- 
nitude lower than typical type lib supernovae and can 
easily be hidden in the spectra. This supernova would 
produce a Ib/c supernova. 

1 1 .35183334 hr 



T500 




-500 



Fig. 11. — SNSPH particle positions for a 5.6 X 10 51 erg explosion 
focused in a 30° wedge of a 25 Mq. The density in this wedge is 
lowered by a factor of 100 with respect to the rest of the star. Shad- 
ing shows the different chemical abundances with black denoting 
hydrogen and lighter shades showing helium, carbon, oxygen and 
silicon abundances. 

How much hydrogen is swept up (and whether or not 
we'll see it in the explosion) depends upon the details 
of the density profile set by the merger ejecta and the 
asymmetry in the explosion. For GRB-like jets expected 
from massive helium-core mergers, the amount swept up 
should not be detectable. But for main sequence stars, 
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the explosion may not be so jet-like, and the amount of 
hydrogen in the explosion could be large. 

Finally, like with most collapsar progenitors, the he- 
lium in the ejecta must either be burned to C/O in the 
shock, or must be hidden in the supernova spectra. This 
remains one of the primary open problems wit h the Col- 
lapsar GRB model as well (jFrver et alj|2007a|) . 

4. COMPARISON TO TRANSIENT POPULATIONS 

We have shown that the rates of our massive he-core 
explosions are on par with the long-duration GRB rate, 
and the total he-merger and main-sequence merger rate 
is roughly 1% the supernova rate. For massive (> 4M©) 
helium cores, energies can be as high as any black-hole 
accretion disk model. These models produce a range of 
hydrogen shells surrounding the star at 1-100 A. U. But 
to truly understand the nature of the explosions of these 
mergers requires a better understanding of the black- hole 
accretion disk engine. In this section, we review some of 
the possible outcomes of these objects. 

4.1. He-Mergers and GRBs 

We have shown that for a subclass of merging systems, 
the accretion rates onto the black hole can provide the 
energetics needed to drive a GRB outburst. In addition, 
the rate of this subclass is sufficient to explain GRBs. 
But can these mergers explain the full suite of GRB di- 
a gnostics. 

iThone et al.1 (|2011l ) argued these mergers could fit a 
very peculiar subset of GRBs with an unusual optical 
counterpart. But can helium mergers explain a broader 
set of GRBs? One of the key constraints of long-duration 
GRB progenitors is the association of a particular sub- 
type of supernova (Type Ic) with GRBs. The 56 Ni esti- 
mates for these supernovae tend to be higher than typi- 
cal Ib/c supernovae with broader line profiles suggest- 
ing asymmetric explosions. Is it possible for the He- 
merger progenitor to reproduce this supernova/GRB as- 
sociation? 

If the helium merger model can not satisfy these re- 
quirments, it may be limited to explaining only a few 
odd GRBs such as the Christmas burst. The shell posi- 
tion in the Christmas Burst (~10A.U.) is characteristic 
of our he-merger predictions, but we expect a range shell 
positions from 1-100 A. U. 

One feature of the helium merger model is that it typi- 
cally has an order of magnitude higher angular momenta 
than the collapsar model. This leads to larger accretion 
disks with larger accretion times, although for the typical 
specific angular momenta of 10 18 cm 2 s _1 , the accretion 
disk remains with 10 9 cm and may not produce drasti- 
cally different characteristics than the Collapsar model. 



At this point, our lack of understanding of the GRB en- 
gine makes it difficult to tie this difference in disk size to 
an observation. But as we better understand the engine, 
this difference may lead to observable consequences. 

4.2. Peculiar Supernovae 

If the explosion is asymmetric but there is too much hy- 
drogen within the supernova explosion to produce a rel- 
ativistic jet, helium helium mergers may produce asym- 
metric Ib/c supernovae or, if they sweep up more hydro- 
gen, some fraction of the type lib supernova population. 
These supernovae would be characterized by large asym- 
metries and, perhaps, larger 56 Ni yields. 

Alternatively, if the disk explosion model produces 
more symmetric explosions, the helium merger model is 
more likely to produce a small subset of the type IIP su- 
pernovae and unless the 56 Ni yield is different, it may be 
difficult to uncover differences between he-merger out- 
bursts from typical core-collapse explosions. Again, it 
is likely that these bursts exhibit stronger asymmetries 
than normal supernovae and evidence of these asymme- 
tries may be the only other characteristic differentiating 
these explosions from normal type IIP supernovae. 

Supernovae have been observed that may exhibit a 
recent shell ejection caused by a helium shell ejection 
(SN2011ht: Romin g et al. 2012). Using the same analy- 
sis as section f3 . 1 . 11 we calculate the outer shell radii for 
both our low-mass helium cores and hydrogen stars 1101 
Note that the distribution has many more shells with 
outer radii below 1-2 A. U. The explosions from these 
mergers may, at peak, look very similar to comparable 
giant stars. A fraction (6% of hydrogen stars, 34% of low- 
mass helium cores) have outer shell radii beyond 10 A.U. 
that may produce the peculiar features of supernovae 
such as SN 2011ht. The broad range of shell separations 
from these systems is in large part due to the fact that a 
sizable fraction of these systems (30% of hydrogen stars, 
17% of low- mass helium stars) arise from systems where 
the compact remnant is kicked into its companion. 

4.3. Radio Bursts 

If, in a GRB-like explosion, the jet sweeps up too much 
hydrogen to become relativistic, the engine does not pro- 
duce much 56 Ni and supernova explosion is asymmetric 
(preventing it from sweeping up much mass), these ex- 
plosions may produce optical weak, radio loud outbursts. 

This project was funded in part under the auspices of 
the U.S. Dept. of Energy, and supported by its con- 
tract W-7405-ENG-36 to Los Alamos National Labora- 
tory. K.B. acknowledges support from MSHE grant N 
N203 404939. 
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